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Abstract 

In this paper we show how the prediction of CMBR angular power spec- 
tra C'i in non-Gaussian theories is affected by a cosmic covariance problem, 
that is (C;,C;/) correlations impart features on any observed C\ spectrum 
which are absent from the average C l spectrum. Therefore the average 
spectrum is rendered a bad observational prediction, and two new pre- 
diction strategies, better adjusted to these theories, are proposed. In one 
we search for hidden random indices conditional to which the theory is re- 
leased from the correlations. Contact with experiment can then be made in 
the form of the conditional power spectra plus the random index distribu- 
tion. In another approach we apply to the problem a principal component 
analysis. We discuss the effect of correlations on the predictivity of non- 
Gaussian theories. We finish by showing how correlations may be crucial 
in delineating the borderline between predictions made by non- Gaussian 
and Gaussian theories. In fact, in some particular theories, correlations 
may act as powerful non-Gaussianity indicators. 
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1 Introduction 



The CMBR temperature fluctuations are often assumed to constitute a 2-D Gaus- 
sian random field. If this were the case it is known that the angular power 
spectrum C l would fully specify the fluctuations' statistics (see [1] for a review). 
While Gaussianity is probably a good working hypothesis in the context of infla- 
tionary scenarios, it is also accepted that in one way or another the fluctuations 
predicted in topological defect scenarios are non-Gaussian ([2, 3, 4]). The issue 
of non-Gaussianity has so far been approached in the form of Gaussianity tests. 
Topological tests [5], peaks' statistics [1], the 3-point correlation function [6], 
and skewness and kurtosis tests [7, 8] have been proposed. If these tests showed 
the fluctuations to be non-Gaussianity then one would have to do more than to 
measure the C l . A whole set of invariants, components of the n-point correlation 
function (n > 2), would then be required in order to fully specify the fluctuations 
([9]). Something less obvious is that non-Gaussian statistics would also affect 
the connection between theoretical and experimental C\ a relationship behind 
any data-analysis strategy. In this paper we show how this may be the case, 
taking as an example a texture low-/ CMBR model [4] known to display strong 
non-Gaussian behaviour. 

Throughout this paper we will use the notation C\ for the angular power 
spectrum of realizations, and C l for its ensemble average. Whereas C\ is a ran- 
dom variable, C l is a number. Since either C\ or C l are essentially sets of com- 
ponents of the 2-point correlation function, they cannot by themselves reflect 
non-Gaussianity. However, the C\ variances involve the 4 t/l moments of the a l m 
distribution, dependent on the statistics. Moreover the C\ are necessarily depend- 
ent random variables in non- Gaussian theories ([9]). As a result, not only will 
the cosmic variance in the C\ be affected by non-Gaussianity (usually in the form 
of a non-Gaussian variance excess), but also the cosmic variance problem 
becomes a cosmic covariance problem. Cosmic covariance can make the C\ 
vs C l comparison troublesome. An example will be given in Section 2.1 showing 
how Ci correlations may impart features on any observed C\ spectra which av- 
erage out to zero in the C l . Whenever this happens the average C l spectrum is 
a bad prediction for the observed CMBR sky. Conversely, the observed C\ also 
becomes a bad estimator for the C\ a fact already hinted at by the abnormally 
large cosmic variance in the C\ found in some texture models known to be very 
non-Gaussian [4]. 

In this paper we approach the cosmic covariance problem with strategies to 
do away with the C\ correlations. In Section 2.1 we show how the cosmic vari- 
ance excess, seemingly connected with correlations, can often be swept under one 
single variable (the random index, say, j/i). We therefore conditionalize the C\ 
spectrum to this variable, and find that in all the conditionalized sub-ensembles 
the cosmic variance is highly reduced and the correlations disappear. An ana- 
logy with a random tilt Gaussian theory proves to describe realistically what is 
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going on. Hence, instead of using C as a prediction for the Ci, we advocate 
the use of the conditionalized spectra C l (\y\) together with the random index 
distribution F(yi). In Sec. 2.2, on the other hand, we apply to the problem a 
principal component analysis [10]. We define a new set of {C}, rotated from the 
original ones, which diagonalize the covariance matrix. In terms of these, realiz- 
ations C\ and averages C l relate in the usual way. We believe that both methods 
constitute new, more sensible strategies for connecting theory and experiment 
in non-Gaussian models. In Section 3 we digress on the implications of what we 
have said on non- Gaussian cosmic variance to the concept of predictivity. Finally, 
in Section 4, we point out the key role played by C\ correlations when confronting 
non-Gaussian theories among themselves and with Gaussian theories. We stress 
the importance of always considering the \C{\ as a whole, and of comparing pre- 
dictions made by different theories in terms of joint C\ probabilities rather than 
marginal distributions. We devise an approximation scheme with which to com- 
pute the cosmic confusion and preference contours in \C{\ space in the presence 
of correlations. We give examples of strong non-Gaussianity indicators, in the 
form of pockets in \C{\ space, where non-Gaussian theories have high probab- 
ility over their competing Gaussian counterparts. We close the paper with an 
outlook of possible applications of what we have said in a more general setting. 
In particular we suggest that cosmic covariance might dramatically undermine 
traditional methods for predicting the defect Doppler peak structure. 



2 Predicting C\ observations in non-Gaussian the- 
ories 

Topological defect observational predictions are often cast in a language borrowed 
from inflation. CMBR defect simulations typically output a number of skies to 
which a spherical harmonics decomposition is applied so as to obtain 

C\ = J2\ a 'm\ 2 - 

m 

Averaging over independent skies one then obtains the angular power spectrum 
C l =< Ci >. Some simulations (eg. [2]) have qualitatively shown an unusually 
high cosmic variance <r 2 (C/). A quantitative formula for the excess variance (rel- 
ative to Gaussian theories) was also provided by an analytical model for textures 
[4]. In Fig. 1 we have plotted the C l spectrum for a texture model according to 
[4]. Superposed are the ±<r/2 cosmic variance error bars for a Gaussian theory 
with the same power spectrum and for the texture theory. Notice how large the 
excess of cosmic variance is for low /. What is the origin of this excess variance? 
Does it imply that the theory is less predictive? Or have we simply applied to 
the theory an inadequate analysis procedure? 
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Figure 1: The C l spectrum for texture Gaussian spots with p = 0.1 and n = 1. 
The inner (outer) error bars are Gaussian (full) cosmic variance error bars. We 
have superposed on it the fitting n l = 1.2 tilted Gaussian theory and associated 
error bars. 



2.1 The random index strategy 

Consider an hypothetical theory in which the ensemble of all Universes can be 
split into sub-ensembles which are Gaussian theories. However, let the average 
C l spectrum in each sub-ensemble be a random variable from the point of view 
of the overall ensemble. A concrete realization of this idea is a tilted spectrum 
Gaussian theory in which the spectral index n l is a random variable. Then 

< Ci{n l ) >„,= C l (n l ) and a 2 nt (Ci(n 1 )) = C l2 (n 1 )^^, where the subscript n l 
means that the averages are taken within a sub-ensemble of constant n l . If fin 1 ) 
is the spectral index distribution one has for the whole ensemble 

C = J dn l fin l ) < Ciin 1 ) > n ,= J dn l C\n l ) =< C\n l ) > (1) 

but now the cosmic variance is 

a\Ci) = j dn*fin%< Cftf) > n , -<7< 2 ) = + §±|^V')) (2) 

in rough analogy with the texture cosmic variance formula. In such a scenario 
the C l spectrum is an unlikely observation. What everyone sees is one of the 
possible C'(n 8 ) spectra, with Gaussian fluctuations around it. Trying to compare 
the observations with the unphysical C l is responsible for the cosmic variance 
excess. It should be obvious that from the point of view of the whole ensemble 
the Ci cannot be independent, as their individual values give away to some extent 
the sub-ensemble they were drawn from. One can check that indeed, although 

< C,(n 8 )CV(n 8 ) > nl = C'(n i )C , '(n i ), one has 

cov(C,, C v ) = J dn l fin l )iC l in l )C l \n l ) - C l C 1 ') . (3) 
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Figure 2: The function P(yi) for n = 1. The bars represent the 68% confidence 
level interval (1.7,3.4). The peak is at 2.66. 



The suggestion is that C\ excess variance, C\ correlations, and the existence 
of a random spectral index are related incidences. A random tilt theory is of 
course less predictive than a fixed tilt theory. Nevertheless, the extra variance is 
essentially <r 2 (n 8 ), not the sum of all the excess variances in the Cj, as one might 
have naively thought. To make the point clear, consider the extreme example of 
a theory with a deterministic tilt, that is, a theory in which each observer sees a 
spectrum of the form 

r, nH + c 



AC) = ^ (s) 



Now let n l be a random variable. A simple minded calculation of the cosmic 
variance leads to 

which is clearly meaningless, as no randomness is seen by any observer. The 
uncertainty in the predictions of such a theory is just <r 2 (n 8 ). 

We will now show that texture theories behave somewhere in between the 
two examples given above, once one realizes that the sky position of the last 
texture, j/i, acts as a random index. The idea is to conditionalize the C\ to y\ 
and see what the statistics are in the conditionalized sub-ensembles. Following 
the notation and results of [4] it is easy to prove that the average C l spectrum 
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Figure 3: The spectrum C l (\y\) for y-y G (1.5,4) for (n = l,p s = 0.1) Gaussian 
spots. 

conditional to y\ is 

C\\y x ) = C[( yi ) + C[( yi ) = ^-^(w h2 ( yi ) + £'dyN(y)W l ' 2 (yj) . (6) 

This is the sum of the last texture brightness plus another term describing the 
average conditional contribution from all the other textures. The y\ distribution 
function is 

P 1 (y 1 )=N(y 1 )e- M ^ (7) 
with M(yi) = Jq 1 dy N(y). It can be checked that indeed 

C'=<C l (\ yi )>= / dy l P{y l )C l {\y 1 ) = ^—^ dy N(y)W"\y) (8) 

JO 47T JO 

in agreement with [4]. We have plotted P\{y\) in Fig. 2 for an = 1 texture model. 
We have also plotted the 68% confidence level interval. Formula (6) suggests that 
regardless of j/i, each observer sees a near white noise regime (C l ~ const) for 
/ up to l c : j (f v +1 } ~ ^(j/i)- The spectrum then crosses over to scale-invariance 
(or slightly tilted). The cross-over l c is a random variable, so the white noise 
feature is absent from the averaged C\ albeit present in any C l (\y\). The C l (\y\) 
spectra have been plotted in Fig. 3 for (n = l,p s = 0.1) Gaussian spots and 
should be compared with C l in Fig. 1. We see that, although C l =< C l (\y\) >, 
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Figure 4: Log-log plot of the C;(|j/i) spectrum for j/i at the peak (middle spec- 
trum) and borders of 68% confidence level interval (top and bottom). The resid- 
ual cosmic variance is plotted as error bars, and the points superposed on them 
correspond to the cosmic variance of a Gaussian theory with a C l (\y\) spectrum. 



the most likely observed spectrum is shaped rather differently from the average 
C l spectrum, as in fact are any of the observed spectra. A spectrum like the 
average C l is a very improbable observation, a fact which is behind the large 
cosmic variances in the unconditionalized C\. 

The residual cosmic variance within each y\ ensemble can be computed by 
introducing the appropriate modifications in the calculation performed in Sec. 3. 3 
in [4] and is 



cr. 



TX 



(C,(\yi)) 




I y y ^dyN(y)W u \y 



U£'dyN(y)W'**(y))\ 



21 + 1 



Cl(yi)C[( yi 



(9) 



In Fig. 4 we have plotted the Ci(\y\) spectrum with ±<r/2 residual cosmic variance 
error bars for y\ at peak value and borders of 68% confidence level interval. 
Superposed on these error bars are the corresponding error bars for a Gaussian 
theory with a C l (\y\) spectrum. We see that the residual variance is either 
comparable to the Gaussian cosmic variance or is smaller. Hence y\ behaves 
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Figure 5: The Pearson's coefficient cor(C/,CV) for texture Gaussian spots with 
p = 0.1 and n = 1. 

like the random spectral index of the examples above, and the conditionalized 
theory contains a residual variance somewhere in between the Gaussian and the 
deterministic examples given. In the next Section we will also show that the 
conditional correlations are negligible, reinforcing the analogy with a random tilt 
theory. Overall, it seems wise to take C l (\y\) combined with P(j/i), and not C\ 
as the real prediction for a texture theory. 

The example given with textures is probably quite general. Non-Gaussianity 
in defect theories arises whenever a particular CMBR feature is due to a single de- 
fect. Conditionalizing the fluctuations to that defect's variables can significantly 
reduce the variance. Some guesswork is always required, but a good recipe is to 
look at the conditionalized cosmic covariance matrix. If there are no residual cor- 
relations left then there are no more hidden random indices. All our predictions 
are then properly parameterized and their variance has been minimized. 

We should say that it is not a tragedy that the predicted spectra depend 
on a parameter which is only statistically predicted by the theory. After all 
the spectral index n\ related to the shape of the inflationary potential, is an 
entirely free parameter in inflationary scenarios (within a range dependent on 
the particular model). The important thing is to realize that the unpredictivity 
introduced by random indices is the random index variance, not the sum of the 
variance surplus in all the unconditionalized C\. 
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Figure 6: Eigenvector components for texture Gaussian spots with p = 0.1 and 
n = 1. The peaks from left to right correspond to C; in decreasing order of cosmic 
variance. 



2.2 C/ covariance matrix and principal components 

An alternative method to get rid of the correlations is to apply to the problem a 
principal component analysis [10]. The idea is to rotate the {C;} frame into a new 
set {Ci} (the principal components) which diagonalizes the covariance matrix. As 
the C\ are uncorrelated, no features may arise in the observed C\ spectra which 
are not present in the average spectrum < C\ >= C l . This procedure does not 
rely on guesswork, is straightforward and general, but the new {C;} are not very 
intuitive. The procedure bears some formal similarities with Gorski's construction 
[11], but the context and motivation are very different. Here we orthonormalize 
the Ci basis with respect to the cosmic covariance matrix of the underlying theory 
(seen as an inner product). In fact our procedure can be applied to any starting 
basis, including Gorski's basis. 

We will exemplify this procedure with textures. The covariance matrix for 
the Ci in all the texture models considered in [4] can be found as a byproduct of 
the calculation in Sec. 4 of [4]: 



cov(CUV) = ^^fl + -P^j) rdyN(y)W l ' 2 (y)W l '*(y) (10) 
47r V < a z > z 



with / ^ /'. Insight into the correlations can be gained from the Pearson's correl- 
ation coefficient 

which always takes values in the range [—1,1], or the relative covariance matrix 

rcov(C / , C v ) = (12) 
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p=0.1 textures r 
n=1.2 tilted spectrum - 



Figure 7: The C\ spectrum for texture Gaussian spots with p = 0.1 and n = 1. 
The inner (outer) error bars are Gaussian (full) cosmic variance error bars. We 
have also shown the fitting n l = 1.2 tilted theory C\ spectrum. 



which factors out the absolute size of the C\. The cosmic covariance matrix 
for textures depends on n and p S} which should hrst be estimated from tilt and 
normalization measurements at intermediate /. In Fig. 5 we have plotted the off- 
diagonal inferior elements of the correlation matrix cor(C/, CV) for our favoured 
texture model. Note how the correlations decay away from the diagonal, and also 
the abnormally high correlations between the low / multipoles. It is curious to 
note that correlations between neighbouring multipoles do not go to zero for high 
/ as fast as one might expect. All correlations become meaningless for large values 
of n or p s . Following [4] it can be proved that the covariance matrix conditional 
to j/i is simply 

cov(C ; ,CV|2/i) = ^^fl + -Ar?V) rdyN(y)W h2 (y)W l ' s2 (y) . (13) 

47T y <C (1 > J Jyi 

We have checked numerically that for all values of y\ in the 68% confidence 
interval, the value of the Pearson's coefficient between adjoining low / falls below 
0.2 (as opposed to 0.8 for the unconditional C\). 

If n and p s are large enough C\ ~ Cj, but as soon as the correlations become 
signihcant the following principal component structure emerges. Arranging the 
C\ in decreasing order of cosmic variance, the hrst C\ is a linear combination 
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Figure 8: V/ for texture Gaussian spots (p = 0.1, n = 1) in the original C\ frame 
and in the uncorrelated frame C\. 



of the hrst few C\ all multiplied by positive coefficients. This C\ summarises 
the low / normalization, and carries an abnormally large cosmic variance. The 
next few C\ are approximately of the form ' y^-'" 1 , that is, they are derivative 

spectra (ss dC l / dl). As we go up in / we recover the C\ ~ C\ regime. The larger 
the p s the earlier we reach C\ ~ C\. In Fig. 6, as an example, we have plotted 
the C\ coordinates in C\ space for a particular texture model (n = 1, p s = 0.1 
Gaussian spots), and one can see the principal component structure mentioned. 
In Fig. 7 we have plotted the corresponding {C 1 } spectrum with Gaussian and 
full cosmic variance error bars. No correlations exist for {C}, so comparison 
with experiment can now proceed in the usual way. There is still an excess of 
variance relative to Gaussian theories but, with the exception of the hrst C\ the 
relative excess is now much smaller (see also Fig. 8). Again it seems that we are 
sweeping all the variance surplus under a single variable whenever we do away 
with the correlations. 



3 Digression on the concept of predictivity 

Besides their data analysis implications, the results in Section 2 also have something 
to say on the controversial issue of predictivity. Cosmic variance makes theories 
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less predictive, but what this means mathematically is not altogether clear. Pre- 
dictivity is a word colloquially used in serious verbal discussions which seldom 
makes its way into the literature. In Gaussian theories the C\ are all independent 
and extensive (additive and positive). Hence it makes sense to define the variance 
of a spectrum as the sum of the C\ relative variances 

a^d) 2 
,JT. C 12 /-/ • 2/ + 1 ' 1 ' 

thus factoring out of the variance expression the absolute size of the spectrum. In 
non- Gaussian theories the C\ are correlated variables and so this definition loses 
its meaning. Loosely speaking it is obvious that correlations reduce the spectrum 
variance from the sum of the C\ marginal variances, but it is not easy to quantify 
this feeling. We may rotate the C\ into an uncorrelated frame {C 1 }, as in Sec. 2.2, 
but then the principal components are not extensive (for instance they may be 
negative). Therefore taking the relative variances is no longer a sensible way to 
factor out the absolute size of the spectrum. In fact not even in Gaussian theories 
does one have a 2 (Ci) oc C 12 . A possible way out of this problem is to compare 
the variance in the C\ with what this variance would be if the original C\ were 
Gaussian random variables 

From these quantities one can then define the spectrum variance relative to a 
Gaussian spectrum variance as 

1 ^max 

V r =nEVi (16) 

with D = / mM -Ln + l' I n this expression correlations have been duly taken into 
account and the spectrum size has been factored out. The predictivity relative to 
the Gaussian predictivity may then be defined as V r = 1/V r . We have computed 
these quantities for the texture p s = 0.1 model and found V r = 1.55. Also in 
Fig. 8 we have plotted V/ for this model in the original C\ frame and in the 
uncorrelated frame C\. We see that we would have grossly overestimated the 
variance of the spectrum by neglecting the correlations. With the exception of 
the first (normalization) variable Cj, the V/ are much smaller in the new frame. 
It turns out that texture models are less predictive than Gaussian theories but 
not as much as overlooking correlations might have suggested. In general the 
lower the p s the less predictive textures are. However we have here considered 
only C\ spectra. One must bear in mind that especially for low p s texture models 
the best predictions might be in the form of multipole shape factors or inter-/ 
correlators [9]. Gaussian theories, on the contrary, are maximally non predictive 
in this aspect. 
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4 Non-Gaussian signals in {C/} space 



In Section 2 we tackled the problem ol predicting C\ spectra subject to cosmic 
covariance. Here we consider how to confront theories with possibly different 
C\ covariance matrices. Clearly correlations cannot be neglected as they may 
constitute the main difference between one theory and another. It is generally 
impossible to find a single variable capturing most of the gap between the two 
theories. Therefore we devise comparison methods which make use of the joint 
distribution functions of the whole \C{\ spectrum. Once more we explain our 
ideas by applying them to concrete examples. We confront the low-/ texture 
models we have used above with tilted spectra Gaussian theories which fit them 
at / G (25,30). As shown in [4] this fit still leaves a significant suppression of 
power at low / in texture models. We now quantify the strength of this signal. 

It will be useful here to use the concept of cosmic confusion between two theor- 
ies T\ and T 2 in a set Q of measurable quantities. This is a measure of the overlap 
of the distribution functions F\(Q) and F 2 (Q) } and is essentially the percentage 
of the two populations which can be put in a one to one correspondence: 

Cq{T u T 2 ) = JdQ mm(F 1 (Q),F 2 (Q)) . (17) 

Cq(Ti,T 2 ) varies between (measuring Q will act as a crucible between the two 
theories) and 1 (T\ and T 2 are the same theory as far as Q is concerned). As 
a guide to the meaning of C we can express it in terms of n-sigmas. This is 
defined as the separation in units of a between the peaks of two 1-D Gaussians 
with variance a which gives an overlap C. More concretely C = erfc(n/(2v / 2))- 
For a 1, 2, 3, 4 a differentiation one has 0.61, 0.31, 0.13, and 0.04 confusion, 
respectively. Given a set Q of variables the confusion between two theories is 
invariant under non-singular transformations on Q. It can also be proved that 
if we ignore one of the variables in the set and marginalize the distributions F\ 
and F 2 with respect to it we can only increase the cosmic confusion. Hence by 
considering the cosmic confusion between texture theories and their fitting tilted 
Gaussian theories in C\ for all 2 < / < 30 (a subset of all the variables) we 
obtain an upper bound on the confusion between the theories. By adding high-/ 
sections of the C l spectrum or m-structure spectra [9] into our predictions we 
can only decrease the confusion. Also a single low-/ variable may exist which 
captures most of the large angle gap between the theories, but the confusion in 
this variable will never be smaller than the confusion in the C\ for all 2 < / < 30. 

Let us now look at \C{\ spectra as points in a D-dimensional vector space 
(where D = l max — l m in + 1)- By varying the free parameters of texture and 
Gaussian models their average spectra {C 1 } span two 2-D surfaces in this space. 
One can set up a map between these two surfaces so as to maximize the cosmic 
confusion (a procedure approximated by the fit performed at / G (25, 30)). Is the 
confusion left by this identification high enough to render low-/ spectra useless? 
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max(F1,F2) — 




Figure 9: The functions max(F tx (C2 } C3), -Fg(C2, C3)) (top) and 
Pref (F tx (C2, C3), Fg(C2, C3)) (bottom) with (right) and without (left) texture 
correlations switched on. The base space lines are contours of iso-preference in 
the C2-C3 plane. Correlations are essential for bringing out the texture preference 
summit (compare the relative heights of the texture and inflation peaks on the 
left and on the right). Correlations also reconfigure the preference contours and 
turn the valley of indecision (Pref = 1/2) into a gorge. 
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Table 1: Cosmic confusion in various sets of variables between four texture theor- 
ies and their fitting tilted spectra Gaussian theories (fit performed at / £ (25, 30)). 

Plots like Fig. 1 are misleading as they ignore C\ correlations and suggest that 
cosmic variance error bars are hypercubes in the \C{\ space. In fact C\ error 
bars are always hyperovaloids, D — 1 dimensional surfaces of equal probability 
inside which a given percentage of the population lives. The principal axes of the 
ovaloid are parallel to the cartesian axes only when the C\ are uncorrelated, but 
even then the ovaloid axes dimensions are not the marginal variances as plotted in 
Fig. 1. We are therefore dealing with a D-dimensional problem which can never 
be factorized into D one-dimensional problems. Due to the complexity of the 
problem we decide here to truncate the analysis at the level of the second moments 
of the C\ distribution. In this approximation the C\ distribution is approximated 
by a multivariate Gaussian distribution with the theory's covariance matrix. The 
ovaloidal error bars become ellipsoidal and one must extend the C\ range to 
C\ £ ( — 00,00). This is a rough approximation, not true even for Gaussian 
theories, for which the joint C\ distribution is a product of X2/+1 1-D distributions. 
However, this approximation does put texture and Gaussian theories at the same 
level of approximation while allowing for texture non- Gaussian features to be 
included, in the form of off-diagonal elements in the covariance matrix. In Table 
1 we show the cosmic confusion between various texture models and their fitting 
tilted Gaussian models in various sets of C\ for 2 < / < 11, computed in this 
approximation. In all cases the confusion in C\ for 11 < / < 30 is nearly 1. 
Although the largest single C\ differentiation between the two theories occurs for 
C2, one has to take the first ten C\ into account to achieve a 2-sigma signal. 

Once a measurement is made the situation changes. Imagine first that one 
makes a measurement of the set Q with infinite experimental accuracy. If F\ (Q) > 
F2(Q) we can then say that the probability of theory T\ over theory T 2 is 

max(F!,F 2 ) 

Pref 12 (Q) = — , (18) 

r\ + r 2 

the preference function. The confusion after the measurement is now 
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The preference function foliates the {C;} space into preference contours. Prefi 2 = 
1/2 represents the valley of indecision. Summits represent the most conclusive 
possible results. There should be at least two summits, one for each theory. 
Imagine now that a measurement is made giving an intermediate / best fit at 
n l = 1.2 for inflation, p s = 0.1 for textures. Let us examine how the low C\ 
could act as a referee. For graphical purposes we have confined ourselves to the 
C2 — C3 plane, and in Fig. 9 we have plotted the preference function, and its 
iso-contours. Notice how different these are if one switches off the correlations. 
It turns out that correlations are essential for pulling the texture summit up to 
the same height as the inflationary summit. Also the valley of indecision has 
much steeper surrounding slopes if texture correlations are taken into account. 
Since the prominent signal associated with textures is only there because of the 
correlations, this signal is to be seen as a sign of non-Gaussianity. 

Naturally, the outcome of the experiment is a domain of results with an 
experimental distribution function. Therefore given a particular experiment, the 
probability of one theory over another (or the experimental cosmic confusion) is 
obtained by integrating expression (18) weighted by the experimental distribution 
over the whole domain. 

5 The message 

This paper is intended as a warning. By custom and tradition the average C l 
spectrum is always taken as a sensible prediction for what the observed C\ should 
be. We have shown with examples how non-Gaussian statistics may render the 
C l a naive and misleading prediction for the observed sky. If you do not have any 
reason to postulate Gaussianity then we recommend the following recipe. Start by 
computing the C l and then always compute the C\ variance. Compare it with the 
Gaussian variance. If a large excess variance exists, then the standard recipe does 
not work. You are dealing with a cosmic covariance problem, and in this paper 
we gave some alternative recipes for making predictions and confronting theories. 
A topical example is defect Doppler peaks [12]. Suppose that an excess variance 
is found for the relevant C\. Then, even if the average C l shows a single bump, 
this might have little to do with what any observer sees. It could happen that 
each observer sees a rich peak structure, similar to the inflation peaks, but the 
peaks' position and height could be random variables (in fact, random indices). 
That being the case, the average C l would be the statistically weighted envelope 
of all the possible peak curves, and the C l would be blind to the peak structure 
seen by any observer. Computer simulations performed conditional to a value of 
the random indices would show a system of peaks. A calculation averaging over 
the whole ensemble would not. Whatever one's methods, in such a scenario the 
sensible prediction would be the C l spectrum inside each sub-ensemble of fixed 
random indices, the residual cosmic variance, and the random indices distribution. 
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The author feels that this type of behaviour (transposed to the context of sample 
variance) might be the entire point in the controversy surrounding the apparently 
contradictory results given by small angle detections [13]. 
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